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Abstract. For adaptive mixed finite element methods (AMFEM), we first introduce 
the data oscillation to analyze, without the restriction that the inverse of the coefficient 
matrix of the partial differential equations (PDEs) is a piecewise polynomial matrix, ef- 
ficiency of the a posteriori error estimator Presented by Carstensen [Math. Comput., 
1997, 66: 465-476] for Raviart-Thomas, Brezzi-Douglas-Morini, Brezzi-Douglas-Fortin- 
Marini elements. Second, we prove that the sum of the stress variable error in a weighted 
norm and the scaled error estimator is of geometric decay, namely, it reduces with a fixed 
factor between two successive adaptive loops, up to an oscillation of the right-hand side 
term of the PDEs. Finally, with the help of this geometric decay, we show that the stress 
variable error in a weighted norm plus the oscillation of data yields a decay rate in terms of 
the number of degrees of freedom as dictated by the best approximation for this combined 
nonlinear quantity. 
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1 Introduction and main results 

Adaptive methods for the numerical solution of the PDEs are now standard tools 
in science and engineering to achieve better accuracy with minimum degrees of 
freedom. The adaptive procedure consists of loops of the form 

SOLVE ESTIMATE MARK REFINE. (1.1) 
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A posteriori error estimation (ESTIMATE) is an essential ingredient of adaptivity. 
We refer to E [6l 13 III [IH ISS IH [SI HOI for related work on this topic. The 
analysis of convergence and optimality of the above whole algorithm is still in its 
infancy. In recent years, there have been some results for the standard adaptive fi- 
nite element method [39ll43[50l[5ll|2l. In|l27l[33l[l2ll23, convergence analysis 
has been carried out for the AMFEM. 

Let i7 be a bounded polygonal in R^. We consider the following homogeneous 
Dirichlet boundary value problem for a second order elliptic PDE: 

/ -div{AVu) = f inn, 

\ u = on on, ^ ' 

where A E L°°(i7; R^^^) is a symmetric and uniformly positive definite matrix, 
and / G The choice of boundary conditions is made for ease of pre- 

sentation, since similar results are valid for other boundary conditions. In ll23l . 
Carstensen presented an a posteriori error estimate for the mixed finite element 
method of (11.21) . and analyzed its efficiency under some restriction of the coeffi- 
cient matrix. In this paper, we shall prove its efficiency without the restriction, but 
at the expense of introducing data oscillation, which is right a component of the 
new concept of error (the total error) (see ll29ll ). 

To summarize the first main result, let Th, Mh x Lh, {ph,Uh),rih,K,osCh de- 
note the meshes, a pair of finite element spaces, a pair of corresponding discrete 
solutions, the estimators and oscillations in turn. We avoid the assumption that 
A'^^Ph is a polynomial on each element, which is required in [|23l for the proof of 
efficiency of the a posteriori error estimator, and obtain the following efficient esti- 
mates for the Raviart-Thomas, the Brezzi-Douglas-Marini, or the Brezzi-Douglas- 
Fortin-Marini elements 

~ \ \^^^^^iP-Ph)\\l2(n) + Wh^^^^'^iP - Ph)\\l2(^n) + \\u - Uh\\l2(^^) + oscl. 

Secondly, we shall analyze convergence and optimality of the AMFEM of the 
form (11.11) . Here we only concern the stress variable error, which is of interest in 
many applications. 

The convergence analysis of the adaptive finite element method (AFEM) is 
very recent, it started with Dofler [i39il , who introduced a crucial marking, and 
proved the strict energy error reduction of the standard AFEM for the Laplacian 
under the condition that the initial mesh To satisfies a fineness assumption. Morin, 
Nochetto and siebert [|49l ISTll showed that such strict energy error reduction can 
not be expected in general. Introducing the concept of data oscillation and the 
interior node property, they proved convergence of the standard AFEM without 
fineness restriction on To which is valid only for A in (2) being piecewise constant 
on To- Inspired by the work by Chen and Feng [|32ll . Mekchay and Nochetto [fSTIl 
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extended this result to general second elliptic operators and proved that the stan- 
dard AFEM is a contraction for the total error, namely the sum of the energy error 
and oscillation. Recently, Gascon, Kreuzer, Nochetto and Seibert (29] presented 
a new error notion, the so-called quasi-error, namely the sum of the energy error 
and the scaled estimator, and showed without the interior node property for the 
self-adjoint second elliptic problem that the quasi-error is strictly reduced by the 
standard AFEM even though each term may not be. 

However, for convergence of the AMFEM, present woks are done only for the 
Laplacian for the lowest order Raviart-Thomas elements ll27l 1121 and any order 
Raviart-Thomas and Brezzi-Douglas-Marini elements f33l . Since the approxima- 
tion of mixed finite element methods is a saddle point of the corresponding energy, 
there is no orthogonality available, as is one of main difficulties for convergence 
of the AMFEM. In this paper, motivated by the two new notions of the error, by 
establishing a quasi-orthogonality result (see Theorem 14.11) we proved that the 
AMFEM is a contraction with respect to the sum of the stress variable error in a 
weighted norm and the scaled error estimator, which is also called the quasi-error. 

To summarize the second main result, let {7^, (M^, L^), {pk, Uk),rik, osck}k>o 
with divMfc = be the sequence of the meshes, a pair of finite element spaces, 
a pair of corresponding discrete solutions, the estimators and oscillations pro- 
duced by the AMFEM in the A;-th step. We prove in Section 5 that the quasi-error 
uniformly reduces with a fixed rate between two successive meshes, up to an os- 
cillation of data /, namely 

Si, + ivUi < o.\si + ^nl) + cosc2(/, r,), 

where a G (0, 1), 7 > 0, 

£l ■= M'^^^iP-PkMlh^Q) + \ \hkdiv{p - Pk)\\l2^Q), 

and osc(/, Tk) is the oscillation of / over Tk (see Section 2.5). We point out 
here that in some cases, even though the stress variable error is monotone, strict 
error reduction may fail. For instance, when pk = Pk+i and f E Lk, from the 
second equation of (12.31) it follows divpk = — /. Then it holds £| = S^^i- On the 
other hand, the residual estimator r]k := rik{pk, Tk) displays strict reduction when 
Pk = Pk+i but no monotone behavior in general. 

Besides convergence, optimality is another important issue in AFEM which 
was first addressed by Binev, Dahmen, DeVore [13] and further studied by Steven- 
son J57l|, who showed optimality without additional coarsening required in |fT3l. 
Both papers lfT3l [571 are restricted to Laplace operator and rely on suitable mark- 
ing by data oscillation and the interior node property. Gascon, Kreuzer, Nochetto 
and Seibert [1291 succeeded in establishing quasi-optimality of the AFEM without 
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both the assumption of the interior node property and marking by data oscillation 
for the self-adjoint second elliptic operator. 

Since all decisions of the AMFEM in MARK are based on the estimator T]k, a 
decay rate for the true error is closely related to the quality of the estimator, which 
is described by the global lower bound 

Hereafter, following the idea in [|29l . we refer to the square root of right-hand side 
above as the total error [|46l . The lower bound demonstrates that the estimator is 
controlled by the error with weights except up to an oscillation term and one can 
observe the difference between Sk and rfk only when oscillation is large. Further- 
more, from the upper bound < r/^ and osc^, < ril it follows Si + osc^ < rjl. 
This implies that the total error , which is the quantity reduced by the AMFEM, 
is controlled by the estimator. Since the estimator itself is an upper bound for the 
quasi-error, in view of the global lower bound it holds 

£l + oscl^r^l^£l + ^r,l. (1.3) 

In short, the behavior of the AMFEM is intrinsically bonded to the total error, 
which measures the approximability of both the flux p = AVu and data encoded 
in the oscillation term. Note that when A^^ph is a piecewise polynomial vector, 
oscillation will reduce to approximation of the right-hand side term / of (|1.2I) (see 
Section 2.5). In general cases, approximation of data A appeared in osc^ couples 
in nonlinear fashion with the discrete solutions pk- 

In Section 6, we shall introduce two approximation classes and based on 
the total error and oscillation of /, respectively. Using a quasi-monotonicity prop- 
erty of oscillation and a localized discrete upper bound, we prove the following 
quasi-optimal convergence rate for the AMFEM in terms of DOFs by assuming 
the marking parameter 9 E (0,6'*) with < 6'* < 1 (see Theorem 16. 31) : 

(4 + os4)i/2 < c'Q^s, omp, /, A)u + ii/i Ug)(#r^ - #ro)-^ 

The rest of this paper is organized as follows. In Section 2, we shall give 
some preliminaries and details on notations. Some auxiliary results are included 
in Section 3 for later usage. Section 4 is devoted to the analysis of efficiency of the 
a posteriori error estimator The proof of convergence for the AMFEM is placed 
in Section 5. Finally, we shall prove the quasi-optimal convergence rate for the 
AMFEM in Section 6 and give conclusions in Section 7. 
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2 Preliminaries and notations 



2.1 Weak formulation 



By splitting (11.21) into two equations, the mixed formulation is given as 

J — div p = f and p = AVu in 

\ u = on d^l. 



(2.1) 



Since the coefficient matrix A is symmetric and uniformly positive definite, by 
the Lax-Milgram theorem, there exists a unique solution u E Hq(Q) to the prob- 
lem (12.11) . Moreover, the weak formulation of (12.11) reads as: Find {p, u) G 

/7(div,fi) X L'^iVL) such that 



where .^(div,^) := {q e LF'iVtf' : div q E L^(fi)} is endowed with the norm 
given by | |g| := Mlh^n) + ||div q\\l2(^^y and (-, ^o.f^ denotes inner 
product on f2. 

For a given shape-regular triangulation Th of Q into triangles, let and 
denote finite dimensional subspaces of i7(div, il) and -L^(fi), respectively. In the 
step SOLVE a mixed finite element method reads as: Find (ph, Uh) E Mh x Lh 
such that 



where fh is the L^— projection of / over L^. 

It is well-known that existence and uniqueness of the solution of (12.21) hold 
true, and that the discrete problem (|2.3I) has a unique solution when a discrete 
inf-sup-condition is satisfied by the discrete spaces Mh and Lh (cf. (221). So we 
are interested in controlling stress variable error e := p — Ph ^ (div, fl) and 
displacement error e := u — Uh E L'^(^l), and suppose that the module SOLVE 
outputs a pair of discrete solutions over namely, {ph, Uh) = SOLVE{Th)- 

2.2 Mixed finite elements 

We consider some well-known mixed finite elements for the discretization prob- 
lem (12.31) . such as Raviart-Thomas (RT) elements, Brezzi-Douglas-Marini (EDM) 
elements and Brezzi-Douglas-Fortin-Marini (BDFM) elements It52l 12111221 . which 
are briefly described for aU triangle T E Thhy some Di{T) C C(T) and Mi{T) C 
C{TY given in the following table. 



(A V, ?)o,!^ + (div g, u)o,f7 = for all g G //(div, fi), 
(div p, f )o,f7 = -(/, for all v E L^{Vt), 



(2.2) 



{A Ph, qh)o,n + (div qh, Uh)o,n = for all qh E Mh, 
(div Ph, Vh)o,Q = -{fh, Vh)o,n for all Vh E Lh, 



(2.3) 
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Examples for mixed finite elements 



Element 


Mi{T) 


Dl{T) 


RT 




Pi 


BDM 


p'^ 1 


Pi 


BDFM 


{q e p-f_^, : (q ■ u)\gT e Ri(dT)} 


Pi 



Here Pi denotes the set of polynomials of total degree < / and Ri (dT) denotes the 
set of polynomials of degree at most / on each edge of T (not necessary continu- 
ous). 

By using the above sets Mi{T) and Di{T), the discrete spaces Mh and Lh are 
given by 

Mh : ={qhe H{diY, Q) : qtW e Mi{T) for T G %}, 
U : ={vneL\n):Vh\TeDiiT)fovTeTh}. 

2.3 Assumption on Th 

Let 7/i be a shape regular triangulation in the sense of ['345 which satisfies the 
angle condition, namely there exists a constant ci such that for all T G 7^ 

Ci^h'^ < \T\ < ci/i^, 

where hx ■= diam(T), and|T| is the area of T. 

Let Eh denote the set of element edges in Th, Jiv)\E '■= i'>j\T_^)\E — {v\t_)\e 
denote the jump of f G H^iXJTh) over an interior edge E := H T_ of length 
He '■= diam(_E'), shared by the two neighboring (closed) triangles T± G Th, spe- 
cially, J{v)\e := {v\t)\e i-i E = T n d^l . Furthermore, for T G Th, we denote by 
ut the union of all elements in % sharing one edge with T, and define the patch 
of E e Eh by 

OOE := |J{^ eTh-.EcT}. 

Denote Th := IJ^^' ^^'^ • -^^(U'^h) L'^(Xh) be an operator with 

HHUTh) ■■= {v G L^in) : VT G Th,v\T G H\T)}. 

Throughout the paper, the local versions of the differential operators div, V, curl 
are understood in the distribution sense, i.e., in D'(il), namely, div/i,curl/i : 
H\\JThf L\n) and Vh : H\\JTh) L\nf are defined such that, e.g., 
diVhv\T := div(v|T) in D'{T), for all T e %. 

2.4 A posteriori error estimators 

For all E G Sh, let r be the unit tangential vector along E, and {ph, Uh) G Mh x 
Lh be the solution of (12.31) with respect to the triangulation Th- Then the local 



6 



estimator is defined by (see Il24l ) 

Vt,k- = ll^''(/ + divp;,)||22(^) + /i|||curl(A-V/.)lli2(r) 

+ \\h'/^JiA-^Ph ■ r) 1 1^2(9^) + \\h{A-^Ph - VhUh)\\l2^T) 
with < K < 1 and the global estimator is given as 

Vh,n ■= 5^ Vt,k- 
Ten 

Here curl?/' := — for ip = [tpi, 1^2)'^ ■ For convenience we also define the 
stress variable error in weighted norm 

■■= \\A-'/'{p-PH)\\h^n) + \\hdiv{p-p,)\\l,^^y 
Note that in this paper, the Curls of a scalar function are involved as 



Cur 10 := 



dx2 ' dxi 

In ||23]| . reliability of the a posteriori error estimator r^/^ ^ with estimates of the 
stress and displacement variables in weighted norm was obtained under a weak 
regularity assumption on A (see Section 4.2 in [jlSl ). whereas efficiency of rih,K 
was derived by assuming additionally that A~^ph is a piecewise polynomial vec- 
tor. In Section 3, we shall prove the efficiency without this additional assumption 
on the coefficient matrix. But we pay the price to introduce oscillation of data. 

In many applications the stress variable is of interest. We define the local 
estimator for the stress variable error as 

VtM,T): = /i|||/-M|22(y) + /i|||curl(A-V)lli2(r) 

+ hT\\J{A-^Ph-r)\\l,^Q^^, 

and define the global error estimator as 

v\ij>h,Th) ■■= r]\{j)h,T). 

We assume that, for a given triangulation Th and a pair of corresponding discrete 
solutions {ph,Uh) e Mh x L^, the module ESTIMATE for the stress variable 
outputs the indicators 

WTM^T)}Ter, = ESTIMATE{p,,rH). 

Then the estimates of the stress variable error e in a weighted norm are reduced 
to (see [i23J) 

W'^'^Whiu) + l|/^dive||i.(^) < C,r]^rM^TH), (2.4) 
where Ci is a constant independent of the mesh size. In Section 3, we shall show 
efficiency of the estimator rij-^{ph, Th) for the stress variable error in a weighted 
norm. 
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2.5 Oscillation of data 

For an integer n > / + 1, we denote by the L^— best approximation operator 
onto the set of piecewise polynomials of degree < n over T e Th or E e Sh, 
denote hy id the identity operator, and set := id — n^. We define oscillation 
osc^i of data as: 

oT^: = ||/.P>rl(A-V)lli.(n) + l|/^'/'^n+i-/(^-V-r)||i.(r^) 

+ \\hP^{A-^Ph-VhUH)\\h^ny 

For the stress variable error in a weighted norm, convergence and quasi op- 
timality of the AMFEM are involved in the oscillations of the data including the 
right-hand side term /. Then we define the oscillation of data as 

osc\(pH,T) : = /i2||P>rl(A-V,)||i,(^) + /iT||P„\iJ(A-V/.-r)||i.(,^) 

+ WHf - fh)\\l^T) forallTer,. 

Finally, for any subset 7^' C %, we set 

osc\{ph,Th) := osc^^(p/j,T) and osc^ := osc'^^{ph,Th). 

We also define oscillation of / as 

osc^(/,rO :=||M/-A)lli^(n)- 

Remark 2.1. For the estimator of the stress variables, let Tube a triangulation, 
Qh e Mfi be given. By substituting p^ with in the definitions of r]-ji{ph,T) 
and oscrf^{ph,T), we can see that the indicator r)ji{qh,T) controls oscillation 
osc-ji{qh,T), i.e., oscr^{qh,T) < ri-ji{qh,T) for all T G In addition, for the 
stress variables, the definitions of the error indicator and oscillation are fully lo- 
calized to T, which means there holds 777^ {qn, T) = rjj-^ {qn, T) and oscj^ {qn-i T) = 
osC'Ti {qH,T) for any refinement ThofTn with T ^Th^Th cmd qn & Mh- More- 
over, a combination of the monotonicity of local mesh sizes and properties of the 
local —projection yields 

r)Th{(lH,Th) < VTh{(1h,Th) and oscr^{qH,Th) < oscrH(g/f, Tj/) Vg// G Mh- 

We note that in this paper, the triangulation Th means a refinement of Th, all 
notations with respect to the mesh Th are defined similarly. Throughout the rest 
of the paper we use the notation A < B to represent A < CB with a mesh-size 
independent, generic constant C > 0. Moreover, A!=:: B abbreviates A < B < A. 
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2.6 The module MARK 



By relying on Dorfler marking, while only concerning the stress variable error 
in a weighted norm, we select the elements to mark according to the indica- 
tors for the stress variables, namely, given a grid Th with the set of indicators 
{VTh{Ph,T)}t&Th ^iid marking parameter 6 G (0, 1], the module MARK out- 
puts a subset of making elements Ain C Th, i-C-, 

Mh = MARK{{7]rAPH,T)}Teru,rH,e), 

such that Ain satisfies Dorfler property 

VTniPH^MH) >Orirj,{pH,TH)- (2.5) 

2.7 The module REFINE 

In the REFINE step, we suppose that the refinement rule, such as the longest 
edge bisection [|53l[54ll and newest vertex bisection [I56ll47ll48l . is guaranteed to 
produce conforming and shape regular mesh. Given a fixed integer 6 > 1, a mesh 
Th, and a subset A4h C Th of marked elements, a conforming triangulation Th 
is output by 

Th = REFINE{TH,MH), 

where all elements of A^// are at least bisected b times. Note that not only marked 
elements get refined but also additional elements are refined to recovery the con- 
formity of triangulations. Let TZ := TZj-h^Th •= Th/{Th H Th) denote the set of 
refined elements, which means M.h ^Tn^Th- 

In general, the number of these additionally refined elements is not controlled 
by i^Ain, that is to say, ^Th — i^Tn cannot be bounded by CAin with a positive 
constant C, which is independent of Th and may depend on the refinement level. 
On the other hand, by arguing with the entire sequence {7fc}fc>o of refinement, 
Binev, Dahman, and DeVore showed in two dimensions that the cumulative num- 
ber of elements added by insuring conformity does not inflate the total number 
of marked elements [.13il . Stevenson generalized this result to higher dimensions 

m. 

Lemma 2.1. ([58]} (Complexity of REFINE). Assume that % verifies condition 
(b) of section 4 in l^58\l . Let {Tk}k>o be any conforming triangulation sequence 
refined from a shape regular triangulation To, where Tk+i is generated from Tk by 
Tk+i = REFIN E{Tk, Aik) with a subset Aik C Tk- Then there exists a constant 
Cq solely depending on % and b such that 

k-l 

#rfc-#ro <Co5^#7W, forallA;>l. 

j=0 
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2.8 Adaptive algorithm 



We now collect the modules described in the previous sections to obtain the AM- 
FEM of the stress variables. In doing this, we replace the subscript H (or h) 
by an iteration counter called A; > 0. Let To be a shape regular triangulation, 
Vo '■= VTo{Po,To) denote the error indicator onto the initial mesh To, with a right 
hand side / G -L^(fi), a tolerance e, and a parameter 9 E (0, 1]. The basic loop of 
the AMFEM is then given by the following iterations: 

The algorithm for AMFEM 

[TN,{PN,UN)]=AMFEM{ToJ,e,e) 
set fc = 0, rjfe = jjo and iterate 
WHILE r?fe > £ DO 
{l){pk,nk) = SOLVE{Tky, 
(2) {miPk,T)}Ten = ESTIMATE{pk,%); 
('i)Mk= MARKavk {Pk,T)}Ten ^Tk^'^)- 
(4) Tfe+i = REFINE(Tk, Mk); k = k + l. 
END WHILe 
Tn = Tk- 

END AMFEM 

We note that the AMFEM for the stress variables is a standard algorithm in 
which it employs only the error estimator {rj-ji {p^, T) jrerfc , does not use the oscil- 
lation indicators {oscj-^ (pk, T)}TeTk^ does not need the interior node property 
for marked elements. 

3 Analysis of efficiency for estimators 

We devote this section to the analysis of efficiency of the a posteriori error esti- 
mator ?7/i K for the stress and displacement variables in a weighted norm. Herein, 
We avoid the additional assumption that A~^ph is a polynomial vector on each 
element, which is necessary for the proof of the efficiency of rjh^K in E3il for the 
RT, BDM, and BDFM elements. 

Lemma 3.1. Let {ph,Uh) G Mh x Lh be a pair of discrete solutions o/ d2.3P . 
denote the operator defined in Section 2.5. Then, for all T eTh, it holds 

hT\\cm\{A-^Ph)\W(^T) < \\A-'/h\\L2(^T) + /iT||P'curl(A-V)l|L2(T). (3.1) 
Proof. From the triangle inequality, we have 

||curl(A-V)||i.(^) < 2(||P>rl(A-V)||i2(r) + \\Ulcml{A-'p,)\\l^^^). 

(3.2) 
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For all T G Th, let ifj^ denote the bubble function on T with zero boundary values 

1 /2 

on T and < -i/^r < 1, the equivalence of norms | ■ \ \l'^{t) and 1 1 ■ Wl'^ij-) for 
polynomials implies 

||n>ri(A-V.)|li.(^) ^ ||4^'n2curi(A-V.)|li.(^) 

= (^Tn2curl(A-Vh),curl(A-Vh))o,T (3-3) 
-(7/>Tn2curl(A-V/^),P„'curl(A-Vh))o,T. 
From e := p — ph, p := AVu, Stokes theory, and integration by parts, we have 



ipTUlcniliA ^ph) ■ curl(A ^pu) 
= j ^Tn^curl(A"Vft) ■ curl(VM - A'^t) 

= - [ V^Tn^curl(A-V) ■ curl(A-^e) (3-4) 
Jt 

= CurKV'Tn^curKA- V)) " (^"'e) 
< |?/^Tn^curl(v4"Vh) Ui(T) I |v4~^e| 1^2(7). 

A combination of (|3.2I) - (|3.4I) . together with an inverse estimation and the prop- 
erty of L^— projection, yields 

||curl(A-Vh)||i2(^) < (||P>rl(A-V)|U2(^) + /i^i||A-V2e||^2(^)) 

x\\cni\{A~^Ph)\\L2{T). 

This implies the desired result. □ 

Lemma 3.2. Let {ph,Uh) G x be the discrete solutions of A2.3\l . P^+i 
denote the operator defined in Section 2.5. Then, for all E G Eh, it holds 

hf\\J{A-^Pn-T)\y^E) < hf\\P^^,J{A-'p,,.T)\\L^E) 

+ hE\\PnCnT\hiA-^Ph)\\L^u>E) + \ \A~^^^e\\L^{coE)- 

(3.5) 



Proof. For all E G Eh, let ipE denote the bubble function on E with the support set 
ue and < ^Pe < I- Put a := J{A-^ph ■ r). since H^^^ : L'^{E) P„+i(^) is 
an L^— projection operator, where Pn+i{E) is a set of polynomials of total degree 
< n + 1 over E, there exists an extension operator P : C{E) — )■ C{uje) Il59ll60ll 
such that 

Pn^^iaU = n^+icr and W^jEPTl^+MWiu^E) - '^^T\\^l+l(^\\uE)■ (3.6) 
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From the triangle inequality, we obtain 

||0-||i2(£;) < 2(||P^+iO-||^2(^) + ||n^+lO-||i2(£;)). (3.7) 



1 /2 

The equivalence of norms | \ip^ ■ Wl'^^e) and 1 1 ■ | |l2(£;) for polynomials implies 

(3.8) 

From integration by parts, we get 

(V'en^^ia, (t)o,e = / tfjEPUl+ia ■ a = / A" V " curK^'EPn^^.a) 

J E -J I^E 

+ / cni\h{A-^Ph)^EPK+i<^ 

J Up, 



(3.9) 



From Stokes theory, and noticing that A ^ph = Vu — A ^e, we have 

f A-'ph ■ cniliijEPIll+icy) = - [ A-'e-cmliijEPnl+ia). (3.10) 
A combination of (IX9l) and (13.101) yields 

(^i?n^+i^> ^)o,E < \\A-h\\L2^^^)\iJEP'n.l+i(r\HHujE) 

+ \\cuMA-^Ph) WlHue) I H^ePK+i^^WlHue) 
This inequality, together with an inverse estimate and (13.61) . yields 

{iljE'n.l+i(T,a)o,E < {\\A-^e\\L2(^ujE) + hE\\cm\h{A-^ph)\\L^UEh 



(3.12) 



''E W^^n+l'^llL^E)- 

We apply Lemma 13.11 to the above inequality (13.121) and use the property of 
L^— projection to get 

1 — i 

{ipE'nl+itr, (j)o,E < (/ibI |P^curU(A~ V) Wl^^e) + I l^"^e| 1^2(0;^))^^" IWUl^e)- 

(3.13) 

From (13.71) . (13.81) . the property of L^— projection, equivalence of norms Wip^J"^ ■ 
\\l'2{e) and II ■ I |l2(£;) for polynomials, and (13.131) . we obtain 

hT\W\\h(E) ~ {h]l'^\\Pn+l(^\\LHE) + hE\\PnCnYlh{A~^Ph)\\LHcOE) 

The above inequality (|3.14l) implies the desired result (13.51) by canceling one 

||c^||l2(£;). □ 
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Lemma 3.3. Let {ph, Uh) € x be a pair of discrete solutions of M.3^ . 
denote the operator defined in Section 2.5. Then, for all T eTh, it holds 

hT\\A''''ph-VhUh\\L^{T) < hT\\Pn{A'^Ph-yhUh)\\L^{T) 

(3 . 1 5) 



L2(T)- 



Proof. Denote /i := A~^ph — VhUh, and let be the bubble function defined in 
the proof of Lemma eflem 3.1. Since 11^ is the L^— best approximation operator 
onto the set of polynomials of degree < n over T E have 

Mlhi^T) < 2(||P^/i||i2(r)||/i||L2(r) + \\Ulfi\\l2(^T^). (3.16) 

The property of the bubble function ipr indicates 

WKl^WhiT) - W^T^Kl^Wh^T) = i^Tllll^,Kl^)o,T ^^^^^ 
= {ipT^llp, /i)o,T - (^Tn^/i, P2^)o,T- 

Since e = p — ph, e = u — Uh, integration by parts and an inverse estimate lead to 

Jt Jt 

= - [ div(^Tn^/i)e - / A-h ■ ^tKi^ 
Jt Jt 

< l|e||L2(T)|V'rn^/u|Hi(T) + ll^"^e||L2(T)||V'Tn^/i||L2(r) 

< {h^^\\e\\L2[T) + \\A-^^'^e\\L2(T))M\L^T)- 

(3.18) 

A combination of (I3.16I) - (I3.18I) yields 

^T||/i||i2(7.) < (/iT||^n/i|U2(T) + ||e||L2(r) + ||v4~^/^e||L2(2^))||/i||i2(r). (3.19) 
The assertion (13.151) then follows from the above inequality (13.191) . □ 

We now prove efficiency of the estimator ///j ^ by using the above three lem- 
mas. 

Theorem 3.1. Let {p,u) and {ph,Uh) G x be the solutions of M.U and 
^2.3P . respectively, and t]/^ « and osch be defined as in Section 2.4 and 2.5. Then, 
for the estimator of the stress and displacement variables for the RT, BDM, and 
BDFM elements, there exists a constant hidden in <, independent of mesh-size, 
such that 

vIk ^ \ \A'^^^{p-Ph)\\l2^n) + \ \h''div{p-ph)\\l2(^^) + \\u - Uh\\l2^^) + oscl 

(3.20) 
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Proof. Notice that for all T G 7^, it holds 



||/;,''(/ + divp/i)||2,2(T) = ll/i^div e||L2(r) for all < k < 1. 

A combination of Lemmas 3.1-3.3 yields the assertion (13.201) by summing over 
all T G and E e Eh. □ 

Theorem 3.2. Under the assumptions of Theorem 13.71 let ilThiPhjTh), £h, <^nd 
osc/i be defined as in Section 2.4 and 2.5. Then, for the estimator of the stress 
variables for the RT, BDM, and BDFM elements, there exists a constant C2 inde- 
pendent of mesh-size, such that 

C2TiT,XPh,Thf <£l + oscl. (3.21) 

Proof. Since —div e = / + div ph, combining Lemmas I3.1H3.2[ and summing 
over all T E Th and E E Ch, we obtain the desired result (13.211) . □ 



4 Auxiliary results 

In this section, we will give some auxiliary results for convergence and quasi- 
optimality of the AMFEM for the stress variable. 



4.1 Quasi-orthogonality 

Lemma 4.1. Given a function f E Ll{Q), there exists a function q E ifo(n)^ 
such that 

div q = f and \\q\\m(n) < \\f\\L^n)■ 
WQ refer to [|20l l4ll4T]| for detailed proofs of this lemma respectively on smooth 
or convex, non-convex and general Lipschitz domains. 

Lemma 4.2. Let Lh, Mh be respectively the discrete displacement and stress 
spaces given in Section 2.2. Set W := H{div,Q) fl L^iVtY for some g > 2, 
and let Hl^' '^d and _L be respectively the (Q) —projection onto L^, the identity 
operator and L^{Q)— orthogonality. Then there exists an operator Uh : W — )■ 
with the following commuting diagram 

w ^ L\n) 

;n/, ;ni, (4.i) 

such that 
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(I) it holds a local estimate (note that i?^(U 7^)^ fl if (div, f2) C W) 
\\h-\id-Iin)q\\mn) <\q\m(\jTu) for all q e H\[jThf f\H{diY,^l)- (4.2) 

(II) Iih approximates the normal components on element edges with 

Vh{id — Iih)q ■ i^Eds = for all E G e^, Vh € L^, q eW, 

where Ve is the unit normal vector along E. 

For the detailed construction of such interpolation operator Uh and proof of 
these properties, we refer to [i44l l3l l22l . Note that the above commuting diagram 
means 

div{id-Uh)W ± Lh- 

Lemma 4.3. Let Th and Th be two nested triangulations, Ul^ be the L'^{Q) 
—projection onto Lh, and {ph,Uh) G Mh x Lh be the solutions of l \2.3\l . Then 
for any T G Th, there exists a positive constant Cq depending only on the shape 
regularity ofTn, such that 

\\uh - TlL„Uh\WiT) < ^/c'oHT\\A-^'^Ph\W{T). (4.3) 

Proof. Let 11^,^ denote the L^— projection operator over L^- For any T G Th, by 
the definition of L^— projection operator TIl„, we have j^iJlL^^ — TlLjj)uh = 0, 
i.e., iJ^Lf, " ^LH)'^h £ Ll{T). We thus can apply Lemma HU] to find a function 
q G H^iTf such that 

div q = (JIl^ - UL„)uh in T and \\q\\m{T) < - UL„)uh\\L^T)- (4.4) 

We extend q to {{^{iiy by zero, since Uh (or Uh) approximates the normal 
components on elements edge, the second result (//) of Lemma |43] implies that 
(Uh — T^H)q G Mh and supp(n/i — IiH)q ^ T. Noticing that div q G Lh, we have 

(div q - [Ul, - ni^)div q, [Ul, - IIlMo,t = (4.5) 

and 

{uh - {Ul, - nL^)uh, {Ul^ - ni^)div q)o,T = 0. (4.6) 

Since T^Lh^h = Uh, a combination of the first equality of (14.41) . and (I4.5I) - (I4.6I) 
yields 

\\Uh-IlLHUh\\l2(^T) = WO^L^ -'^LnWllhi^T) 

= iO^L,, - ^LH)uh, div g)o,T 

= {uh,{IlL,-UL„)dwq)o,T. (4.7) 
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Using the locality of q, the commuting property (14.11) and (12.31) . we obtain 

{uh, (n^,^ - nL^)div g)o,T = (m/i, (IIl^ - nL^)div g)o,n 

= {uh, div(n/, - nj/)g)o,o 

= -{A-^phA^h-T^H)q\T. (4.8) 
The local approximation (|4.2I) of Lemma 14.21 indicates 

(n^ - nH)9)o,r| < ll^"VllL2(T)(lk - nhg||L2(r) + ||g - nHg||L2(2.)) 

< CaHT\\A~^'^Ph\\L^(T)M\m{Ty (4.9) 

Finally, the desired result (14.31) follows from the second inequality of (|4.4I) and 
dH]) - (EJ). □ 

In order to prove the quasi-orthogonality, we need to introduce a pair of aux- 
iliary solutions. Let fu '■= ^L„f denote the L^— projection of / over Lh, and 
consider the following problem: Find {ph, Uh) € x such that 

{A~^Ph, qh)o,n + (div qh, Uh)o,n = for all qn G Mh, 

(4.10) 

(div ph, Vh)o,n = -Uh, Vh)o,n for all Vh e Lh- 

In fact, the solution (j)h,Uh) of this auxiliary problem may be regarded as 
another approximation to the flux and displacement {p,u). 

Lemma 4.4. Let Th and Th be two nested triangulations, osc{ fh, Th) denote the 
oscillation of fh '■= ^L^f over Th, {Ph, Uh) and {ph, Uh) be the solutions of ( 12.31) 
and ( l4.iOD . respectively. Then there exists a constant Cq depending only on the 
shape regularity ofTu such that 

\\A-^'\Vh-Vh)\\LHn) < v^osc(A,rH.) (4.11) 

Proof. Recall that {p^ — Ph,Uh — Uh) G x Lh satisfies the equations 

iA'^{ph - Ph), qh)o,n + (div qh, Uh - Uh)o,n = for all qh G Mh, 
{div{ph - Ph, Vh)o,n = -{fh - fn, Vh)o,n for all Vh G Lh- 

(4.12) 

According to the above equations (|4.12l) . we choose qh = Ph—Ph and Vh = Uh—Uh 
to obtain 

\\A-^/'^{ph-ph)\\\2^^^) = -{div{ph-ph),Uh-Uh)o,n 

= {fh - ^L„f, Uh - Uh)o,n = {fh - ^L„f, Vh)o,n- 

(4.13) 
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Since C L^, it holds 

for allvH G Lff. This implies 

^Lh'^lJ' = ^Laf and {UlJ - Ul^J, IlL„Vh)o,n = 0. 
From (14.131) . we have 

= J2iinL,-nL„)f,v,-iiL,v,)o,T. ^^■'^^^ 

TeTn 

We apply Lemma 1431 to Vh — Hlm'^Ii in (14.141) and obtain 

WA^Hph ~ Ph)\\l2(^n) ^ ^ Hrllfh- fH\\LHT)\\A'Hph-Ph)\\L^T) 

< VC'oosc{fh,rH)\\A-y\ph-Ph)\\LHn), 

which leads to the desired result (14.1 II) . □ 

We state the property of quasi-orthogonality as follows. 

Theorem 4.1. (Quasi-orthogonality j Given f E and two nested triangula- 

tions Th and Th, let {ph, Uh) and {pn, uh) be the solutions of( \2.3\l with respect to 
Th and Th, respectively. Then it holds 

{A-\p-pf,),pf,-pH)o,n<C',/'\\A-'/\p-Ph)\\LHn)Osc{U,TH). (4.15) 
Furthermore, for any 6i > 0, it holds 

{I - Si)\\A-^/\p - PH)\\h^^) 
< \\A-'/\p - Ph) I - \\A~'"{PH - Ph) Wh^n) + g^osc2(A, Th). 

(4.16) 

In particular, ifosc{fh, Th) = 0, then it holds 

\\A-'/\p-PH)\\h^n) = (4.17) 
Proof. Let {ph, Uh) solve the problem (14.101) . then we have 

{A-^{p-ph),Ph-PH\n = -{dw{ph-PH),u-Uh)o,n 

(4.18) 

= ifn - fH,u - Uh)o,n = 0. 
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From the above identity (14.181) and Lemma l44l we obtain 

iA-^{p-Ph),Ph-PH)o,n = iA-\p-ph),Ph-Ph)o,n 

< \\A-'/'{p - p,) \\A-'/'{pH - Ph) Wl^u) 

< Cl^'\\A-y^ip-p,)\\L^n)Oscih,rH), 

(4.19) 

which implies the first result (14.151) . 
Furthermore, notice that 

\\A-'/'iP - Ph) WUn) = U-"\p - Ph) Wh^n) ' U-"\ph - Ph) Wh^n) 

~2{A-\p-ph),Ph-PH)o,n, (4.20) 

then for any 5i > 0, from (14.191 ) and Young's inequality we have 

\\A-'/\p - ph) Wh^a) < U-'^'iP - ph) WUn) ' U-'^\Pk - Ph) Wh^n) 

+6, \\A-'/\p - p,) I li.(^) + f osc2(/,, Th), 

which implies the estimate (14.161) . 

In particular, if osc(//t, Tff) = 0. then from (14.191) it follows {A^^{p—ph),Ph — 
PH)o,n = 0. This, together with (I4.20L yields the relation (14.171) . □ 

Although the oscillation of fh over the triangulation Th appears in the estimate 
of quasi-orthogonality, it is dominated by osc(/, 7^). We refer to [l33l for the 
proof of the following observation. 

Lemma 4.5. Let fh denote the L"^— projection of f over Lh, then it holds 

osc(/h,Tff) < osc(/, Tff). 

4.2 Estimator and oscillation reduction 

In this subsection, we aim at reduction of estimator and oscillation. To this end, 
we relate the error indicators and oscillation of two nested triangulations to each 
other. The link involves weighted maximum-norms of the inverse matrix, A^^, of 
coefficient matrix A and its oscillation. 

For a nonnegative integer m = n — I, any given triangulation Th, and v E 
L°°{^1), we denote by H'^v the best L°°(r2)— approximation of v in the space of 
piecewise polynomials of degree < m, and denote by ut the union of elements in 
Th sharing a edge with T. We further set 

n>:=0, P^v:=itd-U^)v, 
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os4^(A-\T) :=ff|(||P-_,Curl A-i||ioc(^) + if^^||P^A-i||i.(^^)). 
Noticing that is defined elementwise, for any subset C Th finally set 

VTHi^'^'^H) ■= maxr]rjjiA'^,T), osc7-^(A"\ 7^) := max oscr^(A"\ T). 

Remark 4.1. (Monotonicityj The use of best approximation in in the defini- 
tion ofrjj-^{A^^, Th) and oscj-j^{A~^, Th) implies the following monotonicity: for 
any refinement Th ofTn, it holds 

r]rM~\'T'h) < nT„{A~\TH) and oscrM'^^h) < osct„{A^\Th) ■ 

To avoid any smoothness assumptions on the coefficient matrix of PDEs, we 
need to quote a result about implicit interpolation, whose proof can be found in 



Lemma 4.6. (Implicit interpolation) Let ffi and n be two nonnegative integer, and 
uj be either one or two dimension simplex. For a positive integer l we denote by 
: L'^{uj,W) — )■ Prhioj^W) the operator of best L"^— approximation in u, and 
Pi := id-Ul^. Then for all V E L°^{uj,W), V G Pfi(w,M') andrh > n, it holds 

IIP? WIIl^M < Ili^^^-^^^IU^HllV^IU^H- (4.21) 

Lemma 4.7. Let Th be a triangulation. For allT E Th and any pair of discrete 
functions cthjTh G Mh, there exists a constant Ai > depending only on the 
shape regularity of To, the polynomial degree I + 1, and the eigenvalues of A^^, 
such that 

r]Tn{<yH.T) <7]rArH,T) + K,7]rAA~\T)\\A-^/\aH - TH)\Wi^^), (4-22) 

oscru{(JH,T) < oscr„{TH,T) + Aioscr^(A"\T)||A"^/2(aH - th)\\l2{cot)- 

(4.23) 

Proof. We only prove the second estimate (14.231) . since the first one (14.221) is 
somewhat simpler and can be derived similarly. We denote by L'^{Th) the square 
integrable function spaces on Th '■= {J^h- The jump of the tangential component 
defines a linear mapping J : Mh — )■ L'^^ by Jiqn) = J{A~^qH ■ r) for all 
qn G Mh from Mh into Lp^ . Recalling = id — with being the 
L^— projection, denoting qn := oh — th and using the triangle inequality, we 
have 

OSCr^((T//,r) < OiiCrH{TH,T) + HT\\PnCni\{A-^qH)\\L'2{T) 

(4 24) 

+H]!^\\Pl^,J{A-\H-T)\yi&T). 
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We split the curl term as 

curl(A~^g//) = Curl ■ + A^^ : curl q^, 

where Curl A^^ is a vector whose every component is the curl of the correspond- 
ing column vector of A^^, and curl qn is a matrix whose column vector is the 
Curl of the corresponding vector of qn- Invoking Lemma WM with u = T and 
noticing that polynomial degree of gji/ is / + 1, we infer for the first term that 

\\P'^{cnT\A-'-qH)\\LHT) < ||P„~i_iCurl A-1|Uoo(^)||A-i/V||l2(t). (4.25) 

Since curl qfj is a polynomial of degree < /, applying (14.211) again in conjunc- 
tion with an inverse inequality, we obtain for the second term that 

\\P^{A-^ : curl qH)\\L^T) < \\P^_[A~^\\L^(T)\\cnTl qHWiHT) 

^ \\P^^l^^^\\L°°(T)\(lh\m{T) 

< Hj,^ \\P^_iA~^\\lc^{t) I |L2(r)- 

(4.26) 

We now deal with the jump residual. Let T' E Th share an interior edge 
E with T. We write J^A-^qu ■ t) = {{A^^qH)\T - (A"^g//)|T') ■ t and use the 
linearity ofIl'^_^_i, Lemma l431 with u = E, and the inverse inequality | |l2(£;) < 

— 1/2 1 1 

Hrp I Itel U2(r) to deduce that 

\\P^^,{{A'\h)\t ■ r)\\mE) 

= \\{P^^,iA-'qH\T))-T\\mE) < ||P„\i(A-V|t)||l2(E) 

< \\P^lA'^\T\\L°°{E)\\qH\\L^{E) < Hj,^^'^\\P^iA''^\\Lc.(^T)\\qH\\L2(T)- 

(4.27) 

Since Th is shape-regular, we can replace H!p by Ht, a similar argument leads 

to 

\\P^,^,{iA-\H)\T' ■ r)\\E^E) < H-'/'\\P:^^^^^ (4.28) 
A combination of (14.271) and (14.281) yields 

\\P^^,JiA~'qH-T)\\L^E) 

= \\P^^MA-'qH)\T - {A-'qH)\T') ■ r)\\EHE) 

(4 29) 

< ||P„V((^~V)|r ■ r)ll^.(^) + ||P„V((A-V)|t' ■ r)||i.(^) 
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By summing over all edges of element T, from the above inequality (|4.29k we 
get 

\\Pl^,J{A-\H-T)UidT) < //r'^'ll^n"^^"'llLoo(.,)||A-l/2g^||^,(^^). (4.30) 

Finally, the desired result (14331) follows from KIM - KIB and (l430l) . □ 

The following two corollaries are global forms of the above lemma. 

Corollary 4.1. (Estimator reduction) For a triangulation Th with Mh C Th, 
let Th be a refinement ofTu obtained by Th '■= REFIN E{Th, M.h)- Denote 
Ai := 3A^ with Ai given in Lemma 14771 and A := 1 — 2^^/^ > with b given in 
Section 2. 7. Then it holds 



+ (1 + 5-,')^,^l{A~\%)\\A-y\<yn - 
for all (7h G Mh, ah € Mh and any 5^ > 0. 



(4.31) 



Proof. For T G Th, applying the first estimate (14.221) of Lemma l4.7l with an, (^h G 
Mh and using Young's inequality with parameter ^3 > 0, we derive 

Hl + 6^T^ivhi^-\T)\\A-'/\aH-ah)\\h^^^y 

By summing over all elements T E Th and using the finite overlap of patches 
CUT, the above inequality (|4.32l) indicates 

VTMh,rh) < il + 5s)'nl^iaH,Th) 

+ (1 + 53-i)3A?r/^JA-Sr.)||A-V2(^^ _ ^H)\\l.^ny 

(4.33) 

For a marked element T e Mh, we set Th,T ■= {T' G Th\T' C T}. Since 
an G Mh and A'~^ jumps only across edges of %, we have J{A^^aH • r) = on 
edges of Th,T in the interior of T. Notice that 1 1/ - fh\\L'^(T') < 1 1/ - fnWL^iT'), 
we then obtain 

E v'tMhX) < 2-'/YrA^H,T), (4.34) 

since refinement by bisection implies 

liT' = |T'|i/2 < (2-^|T|)i/2 < for all r G Th,T. 

On the other hand, for an element T E Th \ M.h, Remark efrem 2.1 yields 

Vn{(^H,T) < rir„{cTH,T). 
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Hence, from (14.341) and the above inequality, by summing over all T G 7^ we 
arrive at 

ri\{aH.rH) < 2~'/'^'^^{aH,MH)+vU^H,rH\MH) ^^^^^ 

From (14.331) . (14.351) . and the monotonicity rjr^ {A^^, %) < {A^^-, %) stated 
in Remark HTTl we get the desired result (14.311) . □ 

Corollary 4.2. (Perturbation of oscillation) Let Th be a refinement o/Th, and let 
Ai be the same as in Cowllarv \4.1\ Then for all an € M^, ah G M^, it holds 

+2Aios4^ {A-\ To) V2 _ an) | • 

Proof. Remark [2?T] vields oscj-jj {crH,T) = oscr,^{aH,T) for all T G Th ^Th- 
Hence, by the estimate (14.231) and Young's inequality, we get 

os4^(aH,T) < 2os4^K,T) + 2K\osc\{A-\%)\\A-^'\ah - OH)\\l2^^^y 

(4.36) 

By summing over T G Th^Tk and using the monotonicity property osct^^ \ 
Th) < osct^ (^"^, To) stated in Remark l4~n the inequality (14.361) indicates the de- 
sired assertion. □ 



5 Convergence for the AMFEM 

We shall prove in this section that the so-called quasi-error, i.e., the sum of the 
stress variable error plus the scaled estimator, uniformly reduces with a fixed rate 
on two successive meshes, up to an oscillation term of /. This means the AMFEM 
is a contraction with respect to the quasi-error. To this end, subsequently we 
replace the subscripts H, h respectively with iteration counters k,k + \, and denote 
by 

% := V%iPk,Tk) 
the scaled estimator over the whole mesh Tk. 

Theorem 5.1. (Contraction property) Given 6 G (0, 1], let {Tk, {Mk, Lk); {pk, 

Uk)}k>o be the sequence of meshes, a pair of finite element spaces, and discrete 
solutions produced by the AMFEM. Then there exits constants 7 > 0, < a < 1, 
and C > depending solely on the shape-regularity of%, b, rjj-^f^A-^ ,%)> '^"^ 
marking parameter 6, such that 

S'k+i + ivl+i < o^\£l + ivl) + Cosc\f, Tk). (5.1) 
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Proof. For convenience, we use the notations := p — p^, := Pk — Pk+i, 

rjk{Mk) ■.= r]r,{pk,Mk),VoiA~^) ■=VroiA'\%)- 

For any 62 > 0, by Young's inequality and the mesh-size functions hk+i < hk, 
we have 

I l^fc+idiv efc+i| |^2(Q) = I |/?./fc+idiv efc| |^2(j^) — I |/?.fc+idiv |^2(j^) 

-2(/ife+idiv Cfc+i, hk+idiv Ek) 
< IIMivcfcll^^^^^ - \ \hk+idi\ Ek\\l2(^^) 

+(52||/ifc+idiv efc+i||^2(f^) + 52^osc'^{fk+i,Tk), 

which implies the inequality 

(1 - 52)||/ifc+idiv efc+i|li2(f^) < ll/ifediv efc||^2(!^) - \\hk+i<li-v Ek\\\2(^a) 

+5^^osc2(/fc+i,7^), 

(5.2) 

where fk+i ■= 11^^^^/ is the L^— projection of / over Lk+i. 

We combine the quasi-orthogonality (Theorem l4.1l) and (15.21) . and take 82 = Si 
to obtain 

(1 - 6i)£l, < £l - \\A-'''Ek\\l.^^^ + ^^osc2(/,+i,rfc). (5.3) 
Applying the estimator reduction (Corollary 14.11) to (15.31) . we get for any 7 > 0, 



(5.4) 



+^[l + 5-^{i^l-X^l{Mk)} 
+l{l + 5^')K,r^l{A-')\\A'y^Ek\\h^^^ 
+Qflosc\fk+i,Tk). 

In what follows we choose 

7:=l/((l + 53-')Air]o(^~')) (5-5) 

so as to obtain 

(l-Si)£l^,+^vl^, < (l-S,)£l+j{l+S,){r^l-\r^l{Mk)}+^^osc\fk+i,rk). 
By using the reliable estimation (12.41) of the stress variable error, and invoking 
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Dorfler marking property (12.51) . the above inequality yields for any constant a, 

(l-<5i)£:|+i + 7r?i+i < a\l-5,)£l + {l-a\l-5i))£l 

+7(1 + - ^4{Mk)] + ^^^osc\h+i,Tk) 

< a\l - 6i)£i + {1 - a\l - 6i))Cirjl 

+7(1 + 53)(1 - X9')4 + ^^osc\h+uTk) 

Ol 

< a 1(1 - dijtfc H 2 '^fcl 



We choose a such that (1 - a^{l - 5i))Ci + 7(1 + 5-i){l - XO"^) = a^^, which 
indic 3.tc s 

2 (1 - 6,)C, + ^{6,C,h + (1 + ^3)(1 - Ag^)) 
« (l-5i)Ci+7 
We now choose ^3 and 5i such that 

53 < A^V2(1 - \e^) and 6^ < min{l, (l_L^f)M}. 

Then it follows 

5iCi/7 + (l + (53)(l-A^2) < ^i_xe'')5; + {l + 5,){l-\e^) 

< (1-A0^)(1 + ^) = 1, 

which leads to < 1. Finally we set 7 = 7/(1 — 61). Then the desired result 
(ISTTT) follows from (1531) . □ 

We note that the oscillation osc(/, 7^) of the right-hand side term / measures 
intrinsic information missing in the average process associated with finite ele- 
ments, but fails to detect fine structures of /. When the oscillation term osc(/, Tk) 
is marked, it is easy to show the following convergence result. 



Corollary 5.1. (Convergence result) Under the assumptions ofTheorem oJ] there 
exit constants p G (0, 1), 7 > 0, and C > depending solely on the shape- 
regularity of To, b, ?77[){A-i,7b)' ^he marking parameter 9, such that 
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6 Quasi-optimal convergence rate for the AMFEM 



6.1 Auxiliary results 

In this subsection, we aim at the discrete upper bound, which is one key for the 
proof for the quasi-optimal convergence rate. Simultaneously, we shall prove and 
quote some preliminary results. 

Theorem 6.1. (Discrete upper boundj Let % and Th be two nested conforming 
triangulations, {ph,Uh) G Mh x Lh and {ph,uh) £ Mh x Lh be the discrete 
solutions with respect to the meshes Th and Th, respectively, and Th := {T G 
Th '■ T is not included in Th}- Then there exist constants Ci and Cq depending 
only on the shape regularity ofTn such that 

\\A-'/\ph-PH)\\h^n) < CiVt,{ph,Th) + Coosc\fh,TH), (6.1) 

and 

#Th < #Th - #Th. (6.2) 

Proof. The second inequality, i.e., (16.21) . follows from the definition of Th- To 
prove the first one, we introduce the solution {ph, Uh) G Mh x Lh to the problem 
(14.101) . From (12.31) and (14.101) . we obtain div{ph — Ph) = 0, which implies 

{Ph - Ph) ■ vds = 0, (6.3) 

where v is the outward unit normal vector along dVl. Thus ph — pn satisfies 
the conditions of Theorem 3.1 in [43] on the polygonal domain namely it is 
divergence-free and fulfills (16.31) . As a result, there exists ^h G H^{^1) with 
Ph— Ph = Curl iph- Since Ph— Ph ^ Mh, this leads to 

^h G Si+^ := {^h G cm : Curl ^Jh G Mh, ^JhW G Pi+2{T) for all T G %} 

(the definition of Sj^^ is analogous). From (14.101) with Qh = Ph — Ph, we get 

\\A-^/^iPh-PH)\\l2(^^) = iA-\ph-pH),Ph-PH)o,n 

= iA-^ph,Ph- PH)o,n - iA-^pH,Ph- PH)o,n 

= -iA-^pH,Ph- PH)o.n = -{A-^ph, Curl i!h)o,n 

(6.4) 

Since div(Curl^H) = for any ipn ^ ^^h"^, from (12.31) with qn = Curl ipn, 
we have 

{A-^Ph, Curl i)H)o,n = -(div Curl ipn, UH)o,n = for all ipn G (6.5) 
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(6.8) 



To connect Sj^"^ with S*]^^, we need to use some local quasi-interpolation, e.g. the 
Scott-Zhang interpolation, (SS In : Sj^"^ S^j^"^ with 

1 /2 

\\il)h-XHi)h\\L^{E)^Hj^ \il)h\m{LOE) for all E e Eh (6.6) 

and 

W^h - In^hllL^T) < HT\^h\m{L,T) for all T eTh- (6.7) 
We note the quasi-interpolation Zh is local in the sense that if T E Th^iTH or E E 
Ehnen (i.e., T or E is not refined), then {iph-lH^h)\T = Oor {^lJh-THilJh)\E = 0. 

With J'h defined in Theorem 16.11 and ipn = '^H'lph, by using integration by 
parts, a combination of (|6.4I) and (|6.5I) yields 

-p/f)|li2(j^) = -{A-^pH,cnTl{iph-^H))o,n 

= ^ ~ ^'^Ph ■ cmliiph - iPh) 
TeTH 

Applying (16.61) and (16.71) to the identity (16.81) . we arrive at 

-Pi/)|li2(f^) < VTH{VH,J^HMh\H\n) 

< Cl/'vrAPH,J'H)\\A-'/'{p,-pH)\\mn), 

which yields 

\\A-'/\pH-pH)\\LHn) < Cl'STn{PH,:FH). (6.9) 

On the other hand, from Ph — Ph ^ we have 

{A~^{ph-ph),Ph -PH)o,n = -{div{ph -pH),Uh - Uh)o,n = 0. (6.10) 
Then a combination of (16.101) . (16.91 ) and Lemma l44l vields 

\\A-^/^iPh-PH)\\l2(^n) = iA'\Ph-PH),Ph-PH)o,n 

= {A~^{ph -Ph+Ph- PH),Ph -Ph+Ph- PH)o,n 

= \\A~'/\p, - pn) I li.(^) + \\A-^/\pH - ph) Wh^n) 

< Cir]^TAPH,J'H) + CoOsc\U,rH), 
namely the result (16. 1|) holds. □ 
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Remark 6.1. One can also include the second term, osc^ {fh^Tu), of the right- 
hand side in ^(5.71) in the first term rj^^^^pn, J^h)- In doing so, however, one cannot 
expect any relaxation of complexity of analysis, because the oscillation of fh over 
Th still appears in the contraction property (see Theorem \5.1\) . 

Corollary 6.1. (Discrete upper bound) Under the assumption of Theorem \6.1\ 
there exist constants Ci and Co depending only on the shape regularity of Th 
such that 

\\A~'^{ph-pH)\\l2^n) + \M^^iPh-pH)\\l2^n) < Ci7]^^{pH,TH)+CoOSC^{f,TH). 

(6.11) 

Proof. Since 

\\hdw{ph - Ph)\\12^^^ < \\Hdiv{ph - PH)\\l2(^n) = osc^ifh^Tn), (6.12) 
in view of (16.11) and Lemma |431 we obtain the desired result (|6.1 II) . □ 

Note that the constant Cq in Corollary 16. H is actually the constant Cq appeared 
in Theorem 16. II plus 1. Here, for simplicity we still denote it by Cq. 

In what follows we shall prove a stable result in the continuous level. Let fn 
denote the L^— projection of / over Lh, we consider the following problem: 

1 -div(AV^) = inn 
y u = on oil. 

By the Lax-Milgram lemma, there exists unique solution u E Hq^CI) to the prob- 
lem (l6j3l). 

Lemma 6.1. (Stable result) Given a shape regular triangulation Th ofQ, let u G 
HQiVt) and u G HliVt) are respectively the weak solutions to the problems di.2D 
and d<5.i3D . and let p = AVu and p = AVu denote the continuous flux. Then 
there exists a positive constant Cq depending only on the shape regularity of Th 
and the eigenvalues of A such that 

\\A-'/\p-p}\\LHn) < Cl''o^c{f,TH.) (6.14) 
Proof. From Green's formula we have 

U-''\p~P)\\hin) = W"y{u-u)\\l.^^) 

[u - m)V ■ {AW{u - u)) 
= if - fH,u-u)o^n. (6.15) 
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Let C Lh be a conforming finite element space, and set Ut '■= [J{T' G Th '■ 
T' n T 7^ 0} for T G Th- We consider the Clement or Scott-Zhang interpolation 
operator or other regularized conforming finite element approximation operator 
J : Hq (Q) Vh which satisfies 

H^^Wv- Jv\\L2(T)<\\Vv\\L2ij,^) ioi &\\T eTh^v eHliSl). (6.16) 

Existence of such an operator is guaranteed (cf. |l3TJ[551|24l|25l). Recall that 
fn := ^LhI is L^-projection of / over Lh- Th is means (/ - fn, VH)o,n = for 
all vh G Lh- Denoting v := u — u, from (|6.15l) and (16.161) we arrive at 

\\A-'/\p-p)\\hin) = if-fH,v-Jv)o,n 

< E \ w-fH)\\ 

i2(T)||V'y||L2(r) 

< C'J'osc{f,rH)\\A-y\p-p)\\LHn), 
which implies the desired result (16.141) . □ 

We finish this subsection by quoting a counting conclusion from [l57ll for the 
overlay T := Ti (B % of two conforming triangulations 7i and T2, which shows 
T is the smallest conforming triangulation for the triangulations 7i and T2- 

Lemma 6.2. ( Overlay of meshes ) For two conforming triangulations Ti and T2 
the overlay T := Ti (B % is conforming, and satisfies 

#r < #ri + #r2 - #ro. 

6.2 Quasi-optimal convergence rate 

In this subsection, we shall prove the quasi-optimal convergence rate of the AM- 
FEM for the stress variable error in a weighted norm. To this end, we need to 
introduce two nonlinear approximation classes. Let Vn be the set of all trian- 
gulations T which is refined from % and #T < A^. For a given triangulation 
T, let Mj- and Lj- denote respectively the approximation spaces to the flux and 
displacement, {pq-,uq-) G M7- x Lj- be the approximation to (p, m), and hq- be 
mesh-size functions with respect to the triangulation T. The quantity of the best 
approximation to the total error in Vn is given by 

a{N-pJ,A): = MrevAU-'/\p ^ Pr)\\h^n) 

+ \\hrdiY{p - PT)\\Un) + osc^pr, T)}'/^ 
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and for s > we define the nonlinear approximation class as 



A, := {{p,f,A)\ \{pJ,A)l := sup N'a{N-p,f,A) < 00} . 

N>No=#To 

Moreover , the quantity of the best approximation to the right-hand side term / in 
Vn is described by 

||/|Ug:= sup iV^ inf osc(/,r). 

N>No=#To 'G^iv 

By the nonlinear approximation theory [l37l[38]| . we know that if / G then 
I I/I Ug < 00. Here, we recall a result of Binev, Dahmen, and DeVore fill which 
shows that the approximation of data / can be done in an optimal way. The proof 
of this result can be found in lfT3l[T4]| . 

Lemma 6.3. (Approximation of data f) Given an f e a tolerance e, and 

a shape regular triangulation %, there exists an algorithm 

TH = APPROX{f,To,e) 

such that 

oscif, Th) < e, #rH - #ro < c\\f\\'p-'/\ 

We now prove that the approximation pf, generated by the AMFEM concern- 
ing the stress variable converges to p in a weighted norm with the same rate 
— #70) as the best approximation described by A^ up to a multiplica- 
tive constant. We need to count elements added by handling hanging nodes to 
keep mesh conformity (see Lemma [2TI) . as well as those marked by the estimator 
(the cardinality of Aik)- To this end, we impose more stringent requirements than 
for convergence of the AMFEM. 

Assumption 6.1 (Optimality). we assume the following properties of the AMFEM: 

(a) The marking parameter 6 satisfies 9 G (0, 6**) with 

* l + Ci(l + 2Aios4,(A-i,ro))' 

(b) Procedure MARK selects a set M.k of marked elements with minimal cardi- 
nality; 

(c) The distribution of refinement edges on % satisfies condition (b) of section 4 
in /till/. 

The limit value depends on the ratio (6*2/6*1)^/^ < 1, which quantifies the 
quality of approximation to the stress variable of estimator r]ji (pk, %), as well as 
the oscillation osc7^(y4"^, To) of coefficient matrix of the PDEs over To- 
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The following lemma establishes a link between nonlinear approximation the- 
ory and the AMFEM through the Dorfler marking strategy. Roughly speaking, we 
prove that, if an approximation satisfies a suitable total error reduction from Th 
to Th {Th is a refinement of Th), the error indicators of the coarser solutions must 
satisfy a Dorfler property on the set IZ of refined elements. In other words, the 
total error reduction and Dorfler marking are intimately connected. 

Lemma 6.4. ( Optimality marking) Assume that the marking parameter 9 verifies 
(a) of Assumption 6.1, and that f is a piecewise polynomial of degree < I on Th- 
Let Th be an shape regular triangulation ofVL, {pn, uh) G Mh x Lh be a pair of 
discrete solutions of M.3\i . Set /i:=i(l — |2)>0, and let Th, be any refinement 
ofTn such that a pair of discrete solutions {ph,, Uh,) G Mh, x Lh, satisfies 

^h, + osc^h. iPh,,Th,) < niSjj + osc^^(p/^, Th)}- (6.17) 

Then the setlZ := IZru^Th, satisfies the Dorfler property 

VTh^Ph^TV) > 0VTh{ph,Th)- 

Proof. Since / is a piecewise polynomial of degree < / on 7^ , it holds 

osc(/,7^) = and osc{fh,,Th,) = 0, 

which imply 

\\Hdiv{p - pH)\\mn) = and \\Kdiv{p - ph,)\\L2{n) = 0. 

These two relations, together with the lower bound (I3.21h . the condition (16.171) 
and the quasi-orthogonality (14.171) . yield 

{l-2fi)C2v'^^{pH,TH) < {l-2fl){Sjj + OSC?,^{pH,TH)} 

< \\A-y\p-PH)\\h^n) - U-'^'iP-PkMhin) 
+ \\HdiY{p - ph)\\12(^^) - \ \Kdiv{p - ph,)\\l2^^) 
+osc^^ {p H, Th) -2osc\ iPh,,Th,)- 

(6.18) 

We estimate separately the error and oscillation terms. By the quasi-orthogonality 
(14.171) and discrete upper bound (16.11) . we get 



\A-'/Hp - PH)\\h^n) - U-'^Hp - PkMUn) + \\Hdiv{p - PH)\\h^n) 

< c,vUPH,n) 



Kdivip ~ Ph,)\\l2^n) = P ^^"^(Ph, -PH)\\l2(^a) 



(6.19) 



30 



For the oscillation term we argue according to whether an element T E Th 
belongs to the set of refined elements TZ or not. For T E TZwe use the dominance 
osc^j-^{pH, T) < r]^^{pH, T) (see Remark 0)- For T e 7}/ n %„ Corollary gj 
(Perturbation of oscillation) together with oh = Ph and ah = Ph, yields 



< 2A^oscliA-\ro)\\A-'/HPk.-PH)\\l.^ny 



(6.20) 



Combining (16.191 ) and (16.201) we infer that 

osc^^ (ph, Th) - 2osc^^_^ {ph, , Th, ) 

= osc^^{ph, TV) + osc^^(p//, Tn^Th,) 

-2osc\^{ph„TH nTh^) - 2osc\^{ph^,rK\TH) (6.21) 

<osc?^^{pH,n) + 2A,oscl{A-\ro)\\A-^/\ph.-pH)\\h^n) 
< (1 + 2ClAlOs4,(A-^ro))r7^^(pH,7^). 

From (16.181) . (16.191) . and (I6.21L we finally deduce that 

^^-^^-'^^ > 1 + + 2A,T<(-4-^ To)) ^^-^^"' '^"^ = ^^^^«(PH,rH). 

In light of the definitions of 9^, 9, and /i, this concludes the proof. □ 

The fact that procedure MARK selects the set of marked elements Aik with 
minimal cardinality, establishes a link between the best mesh and triangulations 
generated by AMFEM, and forms crucial idea of AFEM (see f5T|). In what fol- 
lows we shall use this fact. 

Lemma 6.5. (Cardinality of Aik) Assume that the marking parameter 9 verifies 
(a) of Assumption 6.1, and procedure MARK satisfies (b) of Assumption 6.1, and 
that f is the piecewise polynomial of degree < I onto %. Let {p, u) solve the prob- 
lem ( 12.71) . and let {Tk, (M^, Lk); {pk, Uk); £k}k>o be the sequence of meshes, finite 
element spaces, the discrete solution produced by the AMFEM, and the stress 
variable error in weighted norm 
If {p., /, A) G As, then the following estimate is valid: 

i^Mk < (l-|)-i/2l(p,/,A)|y^Cy'^{£:| + os4^(p,,rfc)}-V2^ (6.22) 

Furthermore, ifTk+i refinement ofTk obtained by the algorithm REFINE 
with 9 G (0, ^*), then 

^Tk+i - #r. < (1 - A)|y^cf + osc^^(p.,r.)}-v2^ 

(6.23) 
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Proof. Let [e | (p, /, A) \ denote the integer component of e ^/'^ | (p, /, A) | . 

We set e"^ := 4"^C^V(^fc + osc^^(pfc, 7^)), where /i = i(l - |^) > and 

Ca = max{l + 2Klosc\{A~\To),2}, and set iV, := /, A)!^'] + 1. 

Recall that 

^{N, + #ro - /, A) := ^ ^ inf + os4(pr, r)}l/^ 

where 6:f:= p-V2(p _ p^)||2^^^^ + ||/irdiv(p - pr)|li2(f^), 
and 

where is the discrete flux approximation to p := AVu with respect to the mesh 
Te, and is the mesh-size function with respect to %. 

Since there exists % G Vn^+#To-i with 7^ 7^ < iV^ + #7o — 1 such that 

{£l + os4^(p„ r,)}i/2 < (1 + i)a{N, + #ro - /, A), 

where e := min{l, e}. This inequality leads to 

< (iv. + #ro-ir{^l + os4^(p„r.)}i/2 

(6.24) 

< (1 + ~e){N, + #ro - 1) V(iV, + #ro - 1; p, /, A) 
<(! + £) supjv>o iV^^T(iV; p, /, A) = (1 + e) | (p, /, A) |,. 

From the above inequality (16.241) . we obtain 

{£l + osc^^(p„ < + < 2e (6.25) 

and 

#7; - #ro < iV, - 1 < A)|y^ (6.26) 

Let % := %(BTk be the overlay of % and 7fc, /i* denote the mesh-size function 
with respect to %, and (p^,, m^,) be a pair of discrete solutions onto %. We shall 
show that there is a reduction with a factor fx of the total error between and 
Pfc. Notice that 7^ is a refinement of %, and since / is the piecewise polynomial 
of degree < / onto %■ Recall that S^^ := ||y4~^/^(p — P*)|||2(q) + ||/i*div(p — 
P*) I li2(Q)5 by the quasi-orthogonality (14.171) and div(p^, — p^) = 0, we get 

SI+oscUp*,%) = \\A-'^\p-Ps)\\Un) - U''/\P. ~Pe)\\h(^n) 
+ \\Kdiv{p - p^)\\l2(^^^ - \ \Kdw{p^ -Pe)\\h{n) + oscf^ (p*, 7;) 
< £l^+osc'^j-^{p^,%). 

(6.27) 
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By the second inequality (14.231) of Lemma 14771 with G M^^, for all 

T E have 



os4^(p„T) < 2os4(pe,T) + 2A2os4^(A-\r)||A-i/2(p^ _pj||2^ 



Summing on %, the monotonicity of the data oscillation (see Remarks 12.11 and 



lL2(f7) 



(6.28) 



nTUi . we get 

< 2osc2-(p„r,) +2A?os4^(A-i,ro)||A-V2(p^ 

< 2os4^(p„ Te) + 2A2os4^(A-i, To) 

A combination of (|6.25l) . (16.271) . and (16.281) yields 

< Ae^Ca = + osc\{j)k, Tk)}. 

Hence, we deduce from optimality marking (Lemma [64]) that the subset TZ := 
IZj^^T* C Tk verifies the Dorfler property (12.51) for 9 < 9^. The fact that pro- 
cedure MARK selects a subset Ai k C Tk with minimal cardinality satisfying the 
same property (12.51) . and (16.261) leads to 

< #7^<#r.-#rfc<#re + #rfc-#ro-#rfc 
= #r.-#ro<|(p,/,A)|y^e-v^ 

which implies the desired result (16.221) . In the third step above, we have used the 
overlay of two meshes (Lemma [6^ . 

The second assertion (|6.23l) follows from #7a:+i — #7a: ^ #A^fc and the first 
resuk (I021) . □ 

Theorem 6.2. Assume that f is a piecewise polynomial of degree < I onto %, 
then the algorithm AMFEM will terminate in finite steps for a given tolerance 
e. Furthermore, set the algorithm AMFEM terminating in the N—th step, and 
denote by Tn the triangulation obtained in the N—th step. Let {p, /, A) G A^, 
and 0(s, 9) := (1 — ^)^^^'^^ i-qI/^ describes the asymptotics of the AMFEM as 
9 — i- 9^,,0 or s — > 0. Then there exists a constant C, depending on data, the 
refinement depth b, and %, but independent of s, such that 

i^TN - i^To < Ce{s,9)\{pJ,A)\y^e-'/\ (6.29) 
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Furthermore, it holds 

{8% + o^cI^{pn,Tn)Y'' < c^e%s,e)\ipj,A)UifrN - #ro)-^ (6.30) 

Proof. By the contraction property (Theorem 15. II) . and the fact that / is a piece- 
wise polynomial of degree < / onto %, there exists a E (0, 1) such that 

4V + 7^rfe+i(Pfc+i''7;+i) < a^i£i + lV%{Pk,Tk)). (6.31) 

The first assertion is a direct consequence of the above inequality (16.311) . 
Since 

Combining the complexity of REFINE (Lemma [2?T]) and the cardinality (16.221) 
of A^fc, we deduce that 

N-l N-1 

Wn - #ro < 5^ < /3 Y^i^k + osc\{pk, Tk)}-'^'', (6.32) 

k=0 k=0 

where/3:=(l-|)-V2^Cy^^|(p,/,A)|y^ 
From the lower bound (13.211) . we infer that 

8l + ^osc\{pk,Tk) < Sl+^r,\{pk,Tk) < {l + ^){Sl + osc\{pk,Tk)]. (6.33) 

On the other hand, the linear rate a = a{6) < 1 of convergence for the quasi 
error implies that for < A; < — 1 

4-1 + Wt.^.{Pn^i, Tn^i) < a'^'^-'-'^iSl + WtM^ 'T;)}. (6.34) 
We combine the above three inequalities (I6.32I) - (I6.34I) to obtain 

N-l 

^ A:=0 

Since a < 1, the geometric series is bounded by the constant sq = 1 / {1 — a^^^). 
By recalling ri^^_^{pN~u Tn-i) ^ £n-i + IVtn-Ap^-^^ Tn-i), we end up with 

- #ro < se/3{l + ^y/'-Wr^ ,{pn-i, Tn-i)}-'^''. (6.35) 

A combination of the above inequality (|6.35l) and the stopping criteria 

^r^(Piv,r^) (6.36) 

34 



yields the second assertion (I6.29I ). 

By raising the second result (16.291) to the s— th power and recording, we obtain 

e < c'Q%s,e)\{pj,A)UifrN - #ro)-^ (6.37) 

Since {E}^ + osc^^(pAr, Tn)Y^'^ ~ "i]Tn^Vn-, Tn), the above inequality (16.371) and 
the stopping criteria (16.361) imply the third assertion (|6.30l) . □ 

Since some results concerning quasi-optimal convergence rate are obtained 
under the assumption that / is a piecewise polynomial of degree < I onto To- By 
inspiration of these results, we consider the following algorithm which separates 
the oscillation reduction of data / and the quasi error: 

Step 1 [Th, h] = APPROX(/, To, e/2), 

Step 2 [rNAPN,UN)]=AMFEM{THjH,e/2,e). 

The advantage of separating data error and discretization error is that in second 
Step 2, oscillation of data / is always zero since the input data fn is piecewise 
polynomial of degree < / over the initial mesh Th for AMFEM. 

We are now at the position to show the quasi-optimal convergence rate. 

Theorem 6.3. (Quasi-optimal convergence ratej For any f E a shape 

regular triangulation % and a tolerance e > 0. Let 6 G (0, 6^), [T/v, (pat, un)] = 
AMFEM(7^, /, £, 9), and osc% := osc'^^{pn, Tn). If{p, /, A) e A, and f E A'q, 
then there exist a positive constant C depending only on data [A, f), refinement 
depth b, and %, but independent of s, such that 

(4 + os4)i/2<c^e^(.,^)(l(p,/,A)|, + ||/iu.)(#r;,-#ro)-^ (6.38) 

Proof. By Lemma 1631 there exists a triangulation Tk such that 

osc(/,rfc) < £, and #r, - #ro < CWfW'ih'^/^ (6.39) 

Let /fc := n^j./ denote the L^— projection of / over L^, and (p, u) be a pair of 
solutions of (12.11) with respect to the right-hand side term fk- By Lemma [611 
(stable result), we obtain 

\\A-^'\p-p)\Wi^^) < Cl'\sc{f,Tu) < C]!\. (6.40) 

By the definition of A^, and the assumption that (p, /, A) E As, it follows 

{pJk,A) E As and 

\{pJk,A)U<\{pJ,A)\s+\\fU^. (6.41) 
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We then use the stopping criteria (16.361) and apply Theorem [Ol to (p, u) to obtain 

(6.42) 

+ osc^^ (pn, Tn) < Vtn (p^' ^ 

and 

#rAr - #% < Ce{s,9)\{pj,,A)\l/^e-'/^. (6.43) 

Notice that ||/iArdiv(p - p)| |L2(n) < \\h kdw{p - p)\\LHn) = osc(/,7^), and 
e{s,6) > 1. A combination of (I09l) - (l643l) yields 

+ os4^(p^, Tn) < 2\\A-^/\p - p)\\l,^^^ + 2|| Wiv(p -p)||i,(^) 
+2| I A-V2(p _ pjy) 1 12^^^^ + 2| |/;,ivdiv(p - p^) 1 1^^^^^ + osc^^{pn, Tn) 
<2Coe' + 2osc\f,%) + 2e'<e' 

(6.44) 

and 

#7^ — #7o = #7Ar — #7fc + #7fc — #7o 

< C0(s,0)|(p,/,,A)|yVV^ + C||/||i//e-V^ (6.45) 

< CQis,empJ,A)U+\\f\\AO'^^e'y'- 
The above inequality (16.451) implies 

e < c'e'{s,emp, f,A)u + ||/iUg)(#rAr - #ro)-^ (6.46) 

The desired result (16.381) then follows from the above two inequalities (16.441) and 
(16461) . □ 



7 Conclusions 

It is time to recall the main results in this paper. Firstly, we have removed the 
restriction of the coefficient matrix of the PDEs which is required in Carstensen' 
work, and have analyzed the efficiency of the a posteriori error estimator obtained 
by Carstensen, for RT, BDM, and BDFM elements . Secondly, For the AMFEM, 
as is customary in practice, the AMFEM marks exclusively according to the error 
estimator and performs a minimal element refinement without the interior node 
property. We have proved that the sum of the stress variable error in a weighted 
norm and the scaled error estimator, reduces with a fixed factor between two suc- 
cessive adaptive loops, up to an oscillation of data /. This geometric decay is 
instrumental to obtain the optimal cardinality of the AMFEM. Finally, we have 
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shown that the stress variable error in a weighted norm plus the oscillation of data 
yields a decay rate in terms of the number of degrees of freedom as dictated by the 
best approximation for this combined nonlinear quantity, namely we have obtain 
the quasi-optimal convergence rate for the AMFEM. 
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